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1. Introduction 


This paper is devoted to the study of diatomic moiecuiar predissociation 
resonances in the Born-Oppenheimer approximation, at energies ciose to 
that of the crossing of the eiectronic ieveis. In such a situation, we aim to 
provide precise estimates both on the real part and on the imaginary part 
(width) of the resonances. As it is well known, they respectively correspond 
to the radiation frequency and to the inverse of the life-time of the molecule. 


In all of the work, the parameter h stands for the square-root of the inverse 
of the (mean-) mass of the nuclei. The Born-Oppenheimer approximation 
permits to reduce the study to that of a semiclassical system of Schrodinger- 
type operators (see, e.g., |KMSWl IMaMel [MaSoj 1. and the size of the system 
depends on the number of electronic levels that are involved. For instance, 
at sufficiently low energies, this system is scalar, and, typically, one can 
apply the numerous results coming from the semiclassical analysis of the 
Schrodinger operator (see, e.e., |Hal |HeSj| IMal IReSil IHiSil |DiSjt |DyZw| IZwj 
and references therein). 


On the contrary, when several electronic levels are involved, only few results 
are available. One may quote |Bal INal IFLNl iGrMal) . where very particular 
potentials are considered, and |K11 IGrMa21 IMaBr] . where the potentials are 
much more general, but the energy considered is lower than that of the 
crossing. Actually, in this last situation the width of the resonances can 
be estimated by a tunnelling effect through a potential barrier, and it is 
exponentially small (in the parameter h). However, according to chemists, 
these widths correspond to such long life-times that it may even surpass the 
age of the universe! For that reason, it seems more reasonable (and, in any 
case, of interest) to consider situations where the widths are not that small. 

As a matter of fact, this is what is expected when the energy considered 
becomes very close to that of the crossing. However, this case is difficult to 
treat in general, because it corresponds to a somehow degenerate situation 
where the characteristic set of the operator has a singularity of the type 

Here, we study a model with one degree of freedom, where such a phenome¬ 
non occurs. Namely, we consider a 2 x 2 matrix system, the diagonal part of 
which consists of semiclassical Schrodinger operators, and the off-diagonal 
part of a lower order differential interaction. We assume that the two po¬ 
tentials cross transversally at the origin, with value 0, and that, at this 
energy level, one of the two potentials admits a well, while the other one is 
non-trapping (see figure 1). 

For such a model, we study the resonances E = E{h) that have a real part 
0(/i^/^) and an imaginary part 0{h). We actually prove their existence, 
and give their asymptotic behaviour as h —)• 0-|-. In particular, we find that 
their widths behave exactly like except possibly for particular values 
of the limit of h~‘^^^E{h), corresponding to positive zeros of some Airy-type 
function, and for which the width becomes o{h^^^). 
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It is important to notice that, for such a kind of system, none of the standard 
methods used in the scalar case can be applied. For instance, the exact WKB 
method (see, e.g., |FuRal IVb ] ) does not work, because of the presence of 4 
phase-functions j = 1,2 (two for each potentials). Indeed, the method 
requires that, for each of them, there exists at least one direction where 
the real part of difference with the three other ones increases, and here it 
cannot be the case. Even the formal semiclassical WKB constructions can 
be performed only at those points where the two potentials do not cross. 
Finally, the recent method proposed by D. Yafaev in [Ya| (and from which 
this work has been mostly inspired) does not seem to be adaptable to our 
case. 


Therefore, instead of trying to generalise the scalar methods, we have cho¬ 
sen to use particular fundamental solutions of the two scalar underlying 
Schrodinger operators, and to take advantage of the fact that everything 
is known on their behaviours (both semiclassical and at infinity) in order 
to solve the system by an iterative procedure. In that way, we can con¬ 
struct two exponentially decaying solutions on one side, and two outgoing 
solutions on the other side, with good estimates on their behaviour up to 
the interaction point where the two potentials cross. This makes possible 
to compute the Wronskian of these four solutions at that point, and obtain 
in this way the condition of quantization that determines the resonances. 
Then, a careful analysis of this condition leads to precise estimates on both 
the real part and the width of these resonances (see Theorem |2.1[). 


Let us mention that our methods still work for problems on L^(M_|_), with 
potentials Vj(r) behaving like Cjjr^ in 0, with Cj > 0 and 0 < a < 2 (see 
Remark 8.8). Moreover, under an additional assumption on IF, perturba¬ 
tions of size can be admitted, too, such as semiclassical pseudodifferential 
operators of order —2, as it occurs after the Feshbach reduction in the Born- 


Oppenheimer reduction process ( [MaMel IKMSWl IMaSoj ): see Remark 2.2 
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In the next section, we give the precise assumptions we work with, and 
we state our main result. Then, in Section 13.11 we construct fundamental 


solutions to the scalar operators both on the real negative half-line and on 
on the positive half-line, and give estimates on them. These fundamental 
solutions are then used, in Section to construct bounded solutions to the 
system in an iterative way, both on the negative half-line and on the positive 
half-line. In Section these solutions and their derivatives are estimated 
more precisely at the crossing point, and, in Section their Wronskian is 
computed and the quantization condition is written. In order to complete 
the proof, an analysis of this quantization condition is performed in Section 
and additional informations are given concerning the resonant states. 
Finally, an appendix (Section]^ contains recalls on the Airy functions, some 
extensions of Yafaev’s constructions for the scalar Schrodinger equation, and 
a list of formulas that may help the reader. 


2. Assumptions and results 
W e consider a system of Schrodinger equations: 

(2.1) Pu^Eu, P 

where stands for — Pj = -|- Vj{x) (j = 1,2), and W = 

W{x,hD,). 

We suppose the following conditions on the potentials Vi(x), V 2 {x) (see Fig¬ 
ure 1) and on the interaction W{x,hDx)'. 

Assumption (Al) Yi(x), V 2 (x) are real-valued analytic functions on M,s 
and extend to holomorphic functions in the complex domain, 

F = {x G C; |Imx| < (5o(IIex)} 

where do > 0 is a constant, and {t) := (1 -|- 

Assumption (A2) For j = 1, 2, Vj admits limits as Rex —>■ Too in F, and 
they satisfy, 

lim Fi(x) > 0; lim V 2 (x) > 0; 

Rea:^ —oo Recc^ —oo 

a:Gr xer 

lim Vi(x) > 0; lim V 2 {x) < 0. 

Rea:^ + oo Rea:—^ + oo 

xer xer 

Assumption (A3) There exists a negative number x* < 0 such that, 

• > 0 and V 2 > 0 on (—oo,x*); 

• W < 0 < V 2 on (x*,0); 

• V 2 < 0 < Vi on (0, -l-oo), 

and one has, 

V;{x*) =: -TO < 0, Y/(0) =: n > 0, V^{0) =: -T 2 < 0. 

Assumption (A4) IT(x, hDx) is a first order differential operator 
W{x,hDx) = ro(x) + ri{x)hDx, 
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where ro(x) and ri{x) are analytic and bounded in F, and ro(0) ^ 0. 


In this situation, in a neighbourhood of the energy 0, the spectrum of P is 
essential only, and the resonances of P can be defined, e.g., as the values E G 
C such that the equation Pu = Eu has a non trivial outgoing solution u, that 
is, a non identically vanishing solution such that, for some 9 > 0 sufficiently 
small, the function x i-A is in L^(M)©L^(M) (see, e.g., |AgCotiReSi] 1. 

Equivalently, the resonances are the eigenvalues of the operator P acting on 
L'^(S.q) © L^(Me), where Mg is a complex distorsion of M that coincides with 
for X >> 1 (see, e.g., |HeMa | ). We denote by Res(P) the set of these 
resonances. 

For E € C small enough, we define the action, 

rxi{E) 


A{E) := / ^E-Vi{t)dt, 

JxUE) 


IxliE) 

where xl{E) (respectively xi{E)) is the unique solution of Vi{x) = E close 
to X* (respectively close to 0), and it is well-known that, in this situation, 
A{E) is an analytic function of E near 0. 

We also fix Co > 0 arbitrarily large, and we plan to study the resonances of 
P lying in the set Vh{Co) given by, 

(2.2) Vh{Co) := [-Coh2/3, - f[0, Cq/i]. 

For h > 0 and fc G Z, we set, 

-2A(0) + {2k + l)'Kh 


(2.3) 


\k{h) := 


2A'(0)/i2/3 


where A.'(0) = 2 ^|Vi( > 0 is the first derivative of A{E) at 0. Then, 

our result is. 

Theorem 2.1. Under Assumptions (Al)-(A4), for h > 0 small enough, one 
has, 

Res (P) n Vh(Co) = {Ek(h); k G Z} D Vh{Co) 
where the Ek{h) ’s are complex numbers that satisfy, 

(2.4) ReE,{h) = Afc(h)h2/3 _ + 0{h^/^), 

27r^ro(0)^ (/ii(Afc(/i))2 + ^2(Afc(/i))^) -/o n, 

(2.5) ImEkih) = - ^ + 0{h^), 

A (Uj 

uniformly as —)• 0+. Here, the function fii and ^2 are dehned by, 

/‘+00 

gi{t) := / Ai {ry^^{ny-t))Ai{-Ty'^{T2y+ t))dy, 

Jo 

f +00 

h 2 {t) ■■= / A[{Ty^^{T 2 y-t))Ai{-Ty'^{Tiy+ t))dy, 

Jo 

where Ai stands for the usual Airy function. 
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Remark 2.1. One can always choose k = k{h) —>■ +oo in such a way that 
^k{h) —7- po as /i —)• 0+, where po G [—Cq^Cq] is fixed arbitrarily. In this 
case, (2.5) gives 


In particular, if {pi{po), p 2 {po)) 7^ (0,0), then the result produces an equiv¬ 
alent of the width of the resonance as /i —?■ 0+, and it is of order Let 
us observe that one has. 


/ + 00 

Ai - t))Ai {T 2 y + t))dy 

-OO 

= (ti + T2)“^/^Ai (-(rf ^ 


(where the last identity comes from an elementary computation involving 

the Fourier transform of Ai and interpreting pi + pi as a convolution of two 

Airy-type functions) and thus, the possible zeros of pi{po)‘^ + p 2 {po)'^ are 

among those of the function i i—)• Ai (—In particular, they 

are necessarily positive. In the case where ti = T 2 , this phenomenon does 

2 -2 

occur exactly at the zeros of A\{—2iT^ ^t), and for these special values of 
Po the result just says that lmEk(h} is o{h^^^). 

Remark 2.2. Under our assumptions, it can be proved that there exists 
an analytic distortion P 2 of P 2 such that, for z G Vh{Co), one has \\(P 2 — 
z)~^\\c{l‘2) = 0{h~^), while the corresponding analytic distortion Pi of Pi 
satisfies ||(Pi — = C>(dis(z, ct(Pi))“^. Using the ellipticity of Pi 

and P 2 , and denoting by W and W* the distorted operators obtained form 
W and W*, we deduce that, if in addition min|z — Xk{h)\ > 6h for some 
5 > 0 constant, and if sup(|ro| + |ri|) is sufficiently small, then, 

\\K{z)\\ := \\h\P 2 - z)-^W*{Pi - z)-^W\\ < 1/2. 

Therefore, observing that the equation [P — z)u = v is equivalent to ui = 
{Pi — z)~^{vi — hWu 2 ) and {1 — K{z))u 2 = {P 2 — z)~^{v 2 — hW*{Pi — z)~^vi), 
we conclude that if the distance between z G PhiCo) and the eigenvalues of 
P (that are the resonances of P) is at least of order h, then one has, 

( 2 . 6 ) (P-z)-^ = 0{h-^). 


Since the resonances of P are separated from each other by a distance of or¬ 
der h, one can apply the standard perturbation theory and conclude that if B 
is a uniformly bounded operator on L^©L^, then the eigenvalues of P + h^B 
in Ph{Co) differ from those of P by 0{h?‘). In particular, in this situation 
our result remains valid if we perturb P by a semiclassical pseudodifferen¬ 
tial operator h?B = h?b{x, hDx) where b is a bounded analytic symbol in 
r X r. This is typically what happens after the Feshbach reduction in the 
Born-Oppenheimer reduction process: see, e.g., |MaMel [KMSWl iMaSoj . 
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3. Fundamental solutions 

In order not to complicate too much the notations, we write the proof in the 
case Ti = r 2 = 1 only, but it is clear that all our treatment can be performed 
with any other positive values of these parameters. At the end of the proof, 
we explain where the changes occur exactly. 

From now on we fix 0 > 0 small enough, and in the sequel we will use the 
following notations : 

II := (-00,0] = M- ; := Fg{[0,+oo)) = Fe{R+), 

where Fg{x) := x + i9f{x), with / G f{x) = x for x large 

enough, f{x) = 0 for x G [0, Xoo] for some Xoo > 0, and / is chosen in such 
a way that, for any x > Xoo, one has, 

fFeix) _ 

(3.1) Im / - V 2 {t)dt > -Ch, 

JXOQ 

for some positive constant C, where the first integral is taken along (ob¬ 
serve that, for E G 2?/i(C'o) and C > 0 sufficiently large, one always has 
Im y/E — V 2 {t)dt = 0{h)). Performing a Taylor expansion of this 
quantity as 0 —)• 0-|-, and using the fact that, for any A: > 1, one has 
t^V 2 ^\t) —)■ 0 as t —>• -|-oo, we see that it is sufficient to choose x^o suf¬ 
ficiently large and that / satisfies, 

for some 5 > 0 constant. By taking / strictly increasing on [xoo,+oo), we 
see that the only possible problem is near Xoo, but there one can take for 
instance f{s) := 


3.1. Fundamental solutions on I^. On II := (—oo,0], and for E G 
Vh{Co) and j = 1, 2, let be the solutions to {Pj — E)u = 0 constructed 
in Appendix 8.2 (in particular, u~^ decays exponentially at —oo, while 
grows exponentially). Then, the Wronskian y^[u~i\ is independent of 
the variable x and satisfies 


(3.2) 


W[u 


'iT’ 


-2 

vr/ii 


(1 -|- 0(hy) (h —> 0). 


For any k >0 integer, we set, 

'■= {u : /£,—)•€ of class ^ sup |u^-^^(x)| < -|-oo}, 

o<i<A: 

equipped with the norm || u So<i<A:®^P/L we define a 

fundamental solution 

Kj,l : C^g{lL)^C!{lL) (j = l,2). 
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of Pj — E on II by setting, for v G 


■ = 


uJl{x) 


(3.3) 


h^W[ulL,Uj,L\ J-c 


Uj^L(.^)v{t) dt 


+ 




h‘^W[ 




ulL{t)v{t)dt. 


Then, satisfies, 

{Pj - E)K,^l = 1, 

and, because of the form of the operator W, an integration by parts shows 
that we also have, 

Kj^lW, K,,lW* : C{!{Il) ^ C^{Il) (j = 1, 2). 

In view of the construction of solutions to the system, we prove. 
Proposition 3.1. As h goes to 0+, one has, 

(3.4) W hK2,LW* Oihh 

(3.5) II h^K^,LWK2,LW* |lr:(cO(/^))= 0{hl). 

Proof. For j = 1, 2, we set, 

ufrix) := max{|n^^(x)|, \hd^ufj^{x)\}] 


^j,L 


hE 


(3.7) 


(3-6) . _ _ , 

Uj{x,t) := ul^{x)u-^^{t)l{t<x} + UjL{x)ulL{t)'^{t>x} = Uj{t,x). 

Thanks to our choice of iFj,L, and doing an integration by parts, we see that 
we have, 

\hKi^LWv{x)\ = 0{h~3) ( j ?7i(x,t)|n(t)|(it +/iC/i(x,0)|n(0)| 

\hK2,LW*v{x)\ =0{h~3) C/2(x,t)|n(t)|(it +/i[/2(a:,0)|n(0)|^ . 

In particular, 

1 7 ° 2 

(3.8) II/iiF 2 ,LkF* 11= 0(^ s) sup / U 2 {x,t)dt + 0{h3) sxvp U 2 {x,t)). 

x&I]^ J —OO xGli^ 

Moreover, using the asymptotics of and hdxU^^ on II, and fixing some 

constant Ci > 0 sufficiently large, we obtain. 


(3.9) 


If X, t < —Cihi , then, 

U2{x,t) = 0{hl) 


g-|Re/,"(y2-E)i/2|/fe 

\{V 2 {x) - E){V 2 {t) - E)\\ 


If t < —Cih3 < X < 0, then, 

1 g-\Ref*(V2-Ef/^\/h 

U 2 {x,t) = 0{he)—— - 

\{V2{t) - E)\i 

If x,t G [—Cih3,0], then U 2 {x,t) = 0(1). 


(3.10) 
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Hence U 2 {x,t) = 0(1) uniformly, and when x < —5 with <5 > 0 constant, 
there exists a positive constant a such that, 
rO /■-' 5/2 

/ U 2 {x,t)dt= / + = 0(/i). 

J —oo J —oo 

2 

On the other hand, if 5 is chosen sufficiently small and x G [—J, —Ci/ia], 
then, there exists a (different) positive constant a such that, 

r-Cihl 


U2{x,t)dt= / U2{x,t)dt + 0{hi) 

J-25 


r25 -a\ti-\x\i\/h 


= 0(/i3|x|-4) 




1 


-dt + 0{h3). 


Setting t = {hs)3 in the integral, we obtain, 
rO , , „ , 


U2{x,t)dt = 0{h3\x\ 4/i3 6) 


oo I 




ds + 0{h3) = 0{h3). 


Finally, when x G [—01/13,0], we have, 

rO 


»-Cih3 


U2{x, t)dt = 


1-5 


U2{x, t)dt + 0(/i3) 


Thus, we have prove, 
(3.11) 


1 p-at^/h 

= 0{h6) ^ -J —dt + 0{h3) = 0{h3). 

JCihi ti 


rO 


sup / U 2 {x,t)dt = 0{h3), 

X'^0 J—oo 


and, by (3.8) (and the fact that U 2 = 0(1)), (3.4) follows. 

Now, let us prove the estimate for Ml := h? Ki^lWK2^lW* . We see on the 
definition of Ki l and on (3.7) that we have, 

\MLv{x)\<Ch 3 / / Ui{x,t)U 2 {t, s)\v{s)\dsdt 

J —00 J —00 
. rO 


(3.12) 


+ 0/13 


Ui{x, t)U 2 {t, 0)|u(0)|c// 


+ 0/i3[/i(x, 0) / U2{0,t)\v{t)\dt 

J —00 

+ 0/itt/i(x, 0 ) 02 ( 0 ,0)|u(0)I. 


Using (3.11) and the fact that Uj = 0(1) uniformly {j = 1,2), we see that 
the last three terms are 0(/i)sup/^ juj (observe that Uj{t,x) = Uj{x,t) for 
all X, t). 


In order to estimate the integral, we use the following properties of Oi: For 
any <5 > 0 small enough, there exists a > 0 constant, such that. 
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• If x,t < X* — S, then, 


Ui{x,t) = 


• If t < X* — 25 and x ^ [x* — 5, 0], then, 

Ui{x,t) = 0{hh-^/’^); 

• If X G [x* — 4(5,0] and t G [—<5, —Ci/ii], then Ui{x,t) = 0{h^t\~^); 

• If X G [x* — 45, 0] and t G [—Cihi, 0] U [x* — 45, —<5], then Ui{x, t) = 

0 ( 1 ). 


Moreover, by the properties of U 2 , we also know that any part of the integral 
corresponding to |t — s| > (5 with <5 > 0 constant is exponentially small. 

We first consider the case x G (— 00 , x* — 25] for some small positive constant 
(5. 

Then, we see that there exists a constant a > 0 such that, 


r ^0 px*—5 px*—5j2 

// Ui{x,t)U 2 {t,s)dtds = 0{h^) dt 

J J—OQ J—OQ J—00 

+ 0{e-^/^) 




Now, when x G [x* — 25, 0], and still denoting by a every new positive 
constant that may appear, we have, 


0 ^0 ^0 

Ui{x,t)U 2 {t, s)dtds = / dt Ui{x,t)U 2 {t, s)ds + 

—00 Jx*—3S Jx*—4S 


and. 



*-45 


Ui{x, t)U 2 {t, s)ds 


(■-5 


dt 


(■-5/2 


‘-35 


c-Ci/ii 


dt 


l-S 


=0{hi 

+ 0{h^ J 

+ 0{hl) J dt 
1 /-o 

+ 0{h6) I „ dt 


-45 




c-Ci/13 g-o /h 


c-Ci/ii 


1-25 

ro 


/ 2 

l-CihB 




ds 


^-a\t\'2/h 


-ds 




f-Cih^ J-S 


J — ds + 0{h3). 


s 4 
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Hence, 

(3.13) 

rO rO 

/ 

Jx*-35 Jx*-iS 


Ui{x,t)U 2 {t, s)ds = 0 {h 2 ) / ^ dt 

JcihTi 


f.2(S 


'Clhi 


—cults —s2 |/h 


1 1 


-ds 


rS jh 


+ 0{hl) [ ^—dt + 0{hr 

Jcihi i4 


For the first term, the change of variables {t, s) i—)> gives, 


r25 


—\/h 


/Cih 3 Jcih^ 


1 1 


-ds = 0(1) 


s' —a\t—s\/h 


-dtds, 


lc,h 

with C 2 ■= and d' ;= (2(i)^/^. Dividing this integral in two parts, 

depending whether t < s or s < t, and first integrating with respect to the 
larger of the two variables, we obtain. 


(3.14) 


Icihi 


dt 


r-25 


1 2 
Ci/l3 


—olts —s2 \/h 


rS' 


oatjh rS' 


—- ds = 0(1) / dt- I 

t2SS Jc2h Jt 


e-^^/^ds 


= 0{h^/^). 

Moreover, a simple change of variable gives. 


S -at2 fh 

dt = 0{h2). 


I • 

JCih^ t4 


Inserting into (3.13), we deduce that, for x G [x* — 25,0], we have, 

rO 


(3.15) 


Ui{x,t)U 2 {t, s)dtds = 0(/i^/^). 


Finally, going back to (3.12), we conclude 


□ 


3.2. Fundamental solutions on On M_|_, let (resp. be the 

solutions to (Pi — E)u = 0 (resp. (P 2 — E)u = 0) constructed in Appendix 
|8.2[ Setting, 


(3.16) 


U: 


in '■= 2^^e*4(la2 


by Proposition 8.5 we have, 


tt: 


± 

2,R 


(x) ~ {l+ 0 {h))^^^^{E-V 2 {x))-^/^e^^dS ^/E-V 2 {t)dt/h ^ 


Then, all these solutions can be extended as holomo rphi c functions in a 
complex sector around M+, and in particular (thanks to (3.1)), O'= 1,2) 
remains uniformly bounded on and exponentially decaying at infinity 
along /^. Moreover, the Wronskians are independent of the 

variable x and satisfy. 


= ^(1 + (h ^ 0). 

TTAIS 


(3.17) 
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We define a fundamental solution of Pj — E on 


(x) 


dt 


(3.18) , 

uJr{x) r+°° 

I J 


+ 


h‘^y^[Ujji,u^jl] Jx 


where v is in the space of bounded continuous functions on and 

the integrals run over 

Then, using the semiclassical asymptotic behaviour of u^ji on one can 
prove exactly as for (3.4)-(3.5) (xqo playing the role of x*) that we have, 

Proposition 3.2. One has, 

(3.19) II hKi^jiW ||^(^o(je))= 0{h3); 

(3.20) II h^K2,RW*Ki^RW ||£(c0(,«))= 0{h3), 
uniformly as h goes to 0+. 


4. Solutions on II and 

In this section, we use the previous construct ions of fundamental solutions 
in order to construct solutions to the system (2.1) in Ir and 

We first consider the interval Ir. We set, 

0 


WyL ■ = 


^2,L • — 


U. 


2,L 


and we look for solutions u := ^ ^ ^ to (2.1) in Ir, that are close to ^ 

or w 2 I ss, h ^ 0. We rewrite (|2.l|) as, 

(Pi — E)ui = —hWu 2 ‘, 
iP 2 -E)u 2 = -hW*ui. 

Assuming that ui is in we set U 2 = —hK 2 ^LW*ui, so that the system 

reduces to, 

(Pi - P)mi = h^WK2,LW*ui, 

and a solution will be given by any ui G C^{Il) such that, 

Ui = + h^KyLWK2,LW*Ui. 


Now, by Proposition 3.1, the operator Mr := h^Ki^iWK 2 ^lW* is 0{h3) 
when acting on By construction, we also know that is in C^{Il). 

Therefore, we can define. 


(4.1) 


Wl,L ■ = 


Y^j>0 

-hK2,LW* 


e ClUlLf, 
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and wi^L is solution to (2.1) in II, with wi^l ^ as h ^ 0+. In a similar 

way, we can define, 


(4.2) W2,l ■= 


-Ej>oMi{hK,,LWul^) 

+ hK2,LW* E,>o MiihK.^LWul^) 


e C^bilLf, 


that is solution to (2.1) in II, with W 2 ,l —)• l as /i —)■ 0+. 

Remark 4.1. The standard WKB method (see, e.g., |FuRai IVo] ) gives us 
asymptotic expansions for uj^ (j = 1,2) inside (x*, 0), of the form, 

(4.3) 

2/ii 




(E- n(i))-v-sin (AEl+nhl + 1) i + 

^ \ k>l 


2hl 


+ Zll^^E - F.W)-''‘cos ( + ^ 1 

I'K 


hi 


%,L ~ 


\ X 


k>l 


where X 2 {E) is the unique solution of V 2 (x) = E close to 0, i'i{x) := 
T,E)VE-V,(t)dt 


3.2 


Now, on a similar construction can be done by using Proposition 

and by starting from the solutions uj^ (j = 1,2) to {Pj — E)uJj^ = 0 on 

Setting Mr := K 2 ^rW*Ki^rW, by Proposition 3.2 we have ||Mij|| = 
2 * 

0{h3), and thus, defining, 

.443 .^1 ^lR + hKr,RWy^,M=^{hK 2 ,RW*u-y 

V -E,>oMUhK2,RWy^) 


e C^bilR?- 


(4.5) 


_ ( ^^^T^'^j>0^R'^2,R ^ ^(-<0(Te\2 

’"""■“I E,>„mKh 


we see that they are both solutions to (^2.l| on and they respectively 


approach ^:= | ) and 


0 


'^2,R ■ — 


0 

h,R 


, as h —y 0+* 


Remark 4.2. Still by the standard WKB method, we have asymptotic 
expansions for uj^ (j = 1,2) inside (0, + 00 ), of the form. 




-i.- 

712/16 


(4.6) 


^2,R ~ 


(Pi(x) - F;)-V4g-U"i(B) 0i(R^di(i ^ y bpk(x)h^}; 
^ k>i 

(E - V2(x))-^/^el ^-2(B) V^-V2(t)dt^^ ^ y b2,k(x)h'^). 


712 /i6 


k>l 


The solutions we have just constructed are not only bounded, and actually 
we have, 
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Proposition 4.1. The solutions Wj^L given by (4.1)-{4.2), and given 


by (j = 1, 2) satisfy, 

wj,LeL^(lL)eL^ilL) ; r ;,-, rgl 2 ( 4 )© l 2 ( 4 ). 




Proof. It is sufficient to prove that, for any > 1, one has wj^s = 
as |x| —)> oo (on Ir or depending if 5 = L or S' = R). But thanks to 
the exponential decay of Uj{x,t) {j = 1,2) as \x — t\ —)> oo, |x| S> 1, we 


immediately see that (3.11) and ( |3.15[ ) remain valid with Uj{x,t) replaced 
by {x)^Uj{x,t){t)~^. As a consequence, the estimates of Proposition [3^ ex¬ 
tend to the operators (x)^hK 2 ^LW*{x)~^ and {x)^hfKi^iWK 2 ,lW*{ x)~^, 
and since {x)^u~^ G C^{Il) {j = 1, 2), the result for wj^l follows. The same 
arguments apply in and the result for wj^r follows, too. □ 


Now, by general theory on system s of ordinary differential equations, we 
know that the space of solutions to (2.1) that are in Ii (resp. in /^) is at 
most two-dimensional. As a consequence, the previous proposition implies, 


Proposition 4.2. E is a resonance if and only if the four solutions wi^l, 
W 2 ,l, wi^R and W 2 ,r o-re linearly dependent. 


5. Estimates at the crossing point 


In this section, we investigate the asymptotic behaviours of Wj^L{x), Wj^R{x), 
and their first derivative at x = 0. 

We first prove, 

Proposition 5.1. Let xq G (x*,0). Then, for x G [xo,0], one has, 

(5.1) = a±ufj^ + b±u^j^, 

with, 

a- = sin 0{h) ; 6_ = 2 cos -|- 0{h)', 

h h 

a+ = ^ cos + 0{h) ; 6+= - sin -|-C>(/i), 

uniformly as /i —)• 0+. 


Proof. Since (u^^ u^j^) is a basis of solutions to (Pi — E)u = 0, we know 

that (5.1) is verified with, 


a_ = 




u 


l,R^ 


W{ul 






W{uf 


Ri 


We compute these Wronskians at some arbitrary point in x G (x*, 0) , by 
using formula (4.3), Proposition 8.2 and (3.17). By (4.3), we have, 


u^Lyx) = 


2hl 


V^(P-Pi(x))4 


■ cos 


A{E) + vi{x) 


h 


+ 0{hi 


r - vr \ 2h6{E -Vi{x))i . 

(u,y M =- — -.m 


SL>±rM_’Li+o(ht). 

h 4 ' 
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Therefore, using Proposition 8.2 and ( |8.12| ), we obtain, 
2/i“5(^^)l fA{E) + iyi{x) TT 


W(u, T,utr>) = -—i-cos 


h 


-^)Bi'(/i-t6) 


+ 


2h 6{E — Vi)'t . f A{E) + i'i{x) TT 




sm 




h 


-(h-H,) 


+ 0(/l3). 


Then, using the asymptotic behaviour of Bi and Bi' at —oo, and ob¬ 
serving that one has the identity |(—(so that (^()2 =[E — 
this gives. 


_l_ ^ 2h 3 ^ f A{E) + ui{x) vr \ fui tt 


=—— cos 


TT 




h 


-7 COS -- -I- T 


+ 


2h 3 ^ /A{E) + vi{x) vr 


sm 


TT 


h 


-I + + 


2h 3 . A{E) 

sm^^ C>(h3), 


TT 


h 


and thus, by (3.17), 


a_ = sin ^ -|- 0{h). 


The same arguments hold for 6_, a+, 6+, and the result follows. 


□ 


From now on, we set. 


(5.2) d := e^dx, 

and we will use this operator in all the Wronskians that will appear (instead 
of the usual derivative), denoting them by W instead of W. 

Proposition 5.2. For j = 1,2 and S = L,R, there exist complex numbers 
aj^S:/^j,S! such that, 


(5.3) 


wi,s(0) = 

dwi,s(0) = 


^ 1 , 5 ( 0 ) +/ 3 i,s^^rs( 0 ) 

ai^su^giO) 


+ 0(h); 


9ui 5 ( 0 ) 4^/3i,5au]^5(0) 
ai, 5 ^ 11 ^ 5 ( 0 ) 


+ 0(h), 


(5.4) 


IC2,s(0) = 

dw2,s(0) = 


02 , 51 ^^ 5 ( 0 ) 
^ 2 , 5 ( 0 ) + ^ 2 ,sut^s( 0 ) 


+ 0(h); 


a2,sdufg(0) 
du-g(O) + /32,sdu+g(0) 


+ 0(h), 


uniformly as h ^ 0+. 
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Proof. We prove (5.3) and (5.4) in the case S = L only (the case S = R 
being similar). We start with j = 1. By (4.1), we have, 

(5.5) wi^L ■= 


with. 


-f "“pL + 

-n,L 

ri,L := hK2,LW*ufL 


+ 0{h), 


In particular, using the expression (3.3) of we find 

i’i,l(0 ) = -ai,Ln^i(0), 

with. 


(5.6) 


-1 


ai.L := 


hW{u2 Li'^2,L) 4-c 


U2,L(.i)iW*u^^^){t)dt. 


In the same way. 


/iifi,LWri,i(0) =/3i,Ln+^(0), 


with, 

(5.7) 


Pl,L — 


1 




/lVV’(u^ 2,) J-OO 

In addition, we can write, 

ri,Lix) = -ai,Lri^^(x) + ri^L{x) 

with, 

I 


ri,Lix) ■■= 


/i>V(<^,M2,l) Jo 


U2,Li^)u2,Li^) - U2^Li^)u2,L( 

X {W*uf^){t)dt. 


In particular ^(0) = 0, and thus r[ ]^{0) = —ai^i(M^^)'(0). We see in the 
same way that '(/i77i,LlI^ri,L)'(0) = /3 pl(m)^^)'(0). 

Moreover, concerning the remainder term in the derivat ive, we observe that 
for any j. S', one has du^g{0) = 0(1) (see, e.g., Remark 8.3). Therefore, for 
any v G C^(Il), we see on the definitions of Ml and 772,L that we have. 


d(hK2,LW*v)(H) = 0(e-^) 

fO 


W2A't)\-H't)\dt + h\{u2r) (o)IXo)l ; 


' —oo 


dMLv(Q) = 0(e 


- 2 \ 


r-0 


|C/ 2 (t, s)|.|u(s)|dtds +/i / |[/ 2 (t, 0)|.|i;(0)|dt 

J —OO 

nO 


where U 2 is defined in (3.6). Then, we observe that \u 2 Lit)\dt = O(e^) 
and, with the same proof as for (3.15), we have \U 2 {t, s)\dtds = 0{e^). 

Using also (3.11), we obtain, 

d(hK 2 lIU*u)(0) = 0(e) sup |u|; 

II 

dMLv(0) = O(e^) sup |i;|. 

Il 
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As a consequence, 


d{hK^^LW*Y.Miul^m = 0{e^); 


j>i 


i>2 


and (5.3) follows. The proof of (5.4) is almost the same, with the only 
difference that the starting function is hKi^iWu^ instead of But 

thanks to the decay properties of L’ that its behaviour is similar 

to that of and (5.4) follows. 


In addition to (5.6)-(5.7), the other constants appearing in (5.3)-(5.4) are, 
(5.8) 


a2,L = 


«1,R = 


Oi2,R = ^ 


~ /o^°° u-R{t){W*u-j^){t)dt 
h'W{u2 R, U 2 J 1 ) 

- lo’^ u^,Rii)iWu2^j^){t)dt 


g ^ I-oo U 2 ,Li'>^){W*r 2 ,L){t)dt 

fo°° '^T,R{t)iWri,R){t)dt^ 


I 3 i,r = 

1^2,R = 


U2,R{t){W*r2,R){t)dt 
h'W{u2 R, U2 J1) 


hW{u^g^, rti jj) 
where we have set, 

r2,L ■= hKi^LWu2 R ; n^R := hK2,RW*u{j^ ; r2,R := hKi^RWu2 R, 
and where, in the case of aj^R and /3j,R, the integrals run over □ 


Now, setting, 
(5.9) 

we also have. 


^+00 

liA{t) := / Ai (y - t)Ai (y + t)dy] 
JQ 

p +00 

Ai(y-t)Bi(y + t)dy, 

Jo 


Proposition 5.3. As h 0+, one has, 
(5.10) 


AE) ^ 


+ 0(8-^), 


aj^L = — 2e7rro(0) ^/iy 4 (Rep) sin ^ Ays (Rep) cos 

(p^(Rep) - fys(Rep)) + 0{e^), {j = 1,2); 

Im/3i,L = 0{h)-, 

Ivn^i^R = 7r^ro(0)^e^ (^^(Rep)^ + pB(Rep)^) + 0{h) 

= 2e^7r^ro(0)^ 


Ai {t — Re p)Ai {s — Re p) 


0<s<i 


X (Ai {t + Rep)Bi (s + Rep) — Ai (s + Rep)Bi {t + Rep)) dsdt + 0{h)-, 
/32,s = 0{e^), {S = L,R). 
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Proof. Let us first study ai^L given in (5.6). Using (3.2) and the exponential 
decay of away from 0, we obtain, 
fO 


— TT 


ai,L = ^ J + ^A(n%,L)(i))dt + C’(e^), 

where 5 > 0 is arbitrarily small. Further, an integration by parts gives, 



— TT t 

0 

ai,L = 


U2,Li*yoit)u 

5 

and since 

Kuf^L 

y{t)uf^{t) is c 

and C > 

1 large 

enough), and 

(5.11) 


/-o 

2e J_s 


U2,Li^)ro{t)u^^L{t)dt + 0{e^). 

We introduce a large extra-parameter A S> 1, and, dividing the integral in 
two parts, we set. 


2 

_ r—Xh^ 

w L 

rO 




-IT f _ 

lW = U- / 2 % L L^t)dt. 

J-XhTi 


We have, 


o^i(A) = 0(£-i) f , = 0(e£-'"“'5, 


hi 

'Xhi y/i 

where c > 0 is a constant, and the estimate is uniform with respect to e > 0 
small enough and A > 1 large enough such that \h?/^ —)• 0. In particular, 
taking, 

(5.12) A > (c“^|lne|)2/3, 

we obtain, 

al^[X) = 0{e^). 


On the other hand, using Propositions |5.1[ |8.2| and 8.5 we obtain, 
—27r , . Al, 


= —(sin-r) 


^ ' h j-xhi. 

rO 


’’o(i)(^i?2) 2Ai(/i 3^i)Ai(/i ii2)dt 


27r A I 

- w(cosu) / 2 

^ “ J-Xhi 


ro{t){C'A2) 2Bi(/i 3^i)Ai(/i 3 ^ 2 )^^ + C>(e^)- 


Making the change of variable y := e ‘^t, and using that, for y G [—A,0], 
we have ro{e^y) = ro(0) + ©(e^A), ^'(e^y) = 1 + 0{£^\), Ai{e-^f,i{e‘^y)) = 
Ai (y-p)+0(A2£2), Ai {e-\2{e^y)) = Ai {y+p)+0{\^e^), Bi {e-^ii{e^y)) = 
Bi (y-p) + 0(AV), we obtain, 

A 

^i,lW = ~ 2evrro(0)(sin —) J Ai (y - y)Ai (y + p)dy 




- 2e7rro(0)(cos —) / Bi (y - y)Ai (y + p)dy + C>(e ) + C>(AV). 
t/ —A 
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Then, using the behaviour of Ai, Ai and Bi at —oo, this leads to, 

~ f A A 

= -2e7rro(0) / Ai (y + p) (sin —)Ai {y - p) + (cos —)Bi {y - p) ] dy 


. h' 

+ + 0(e2) + 0{X^£^), 


h ' 


with c' > 0 constant. Therefore, taking A := ((c^')“^| lne|)^/^ with c = 
min{c, c'}, and using the fact that Imp = 0(e), we obtain the required ap¬ 
proximation of The approximations of 02 ,l, oi,r and 02 ,r are obtained 
in a very similar way (starting from (5.8)), and we omit the proofs. 

Now we study /3 i,r. We have. 




/ u+^is){W*u^ji)(s)ds 


n +00 


+ulR(t)(W*u^ j^){s)ds / U2 j^(s){W*u^ j^)(s)ds 


dt + Oih), 


where the integrals run over Because of the exponential decay of 
away from 0, we immediately see that only a neighbourhood of 0 contributes 
to the integrals (up to 0{e~^^^) with c > 0). Moreover, making an integra¬ 
tion by parts and using the behaviour of all the functions near 0, we see as for 
ai^L t hat, u p to an error 0{e^), W and W* can be replaced by tq. Finally, 
using (3.16) and the expressions of U 2 and given in Propositions 
8.2 and and proceeding as for ai,L, we obtain. 


/3gR = 2i7r^ro(0)^e^ / / Ai{t- p)Ai (s - p)Aout(t + p)Ain{s + p)dsdt 

lo<s<t 


+ 2z7r^ro(0)^e^ 
+ Oih), 


Ai (t - p)Ai (s - p)Ain(t + p)Aout(s + p)dsdt 


0<t<s 


where we have set, 

Aout •— Ai zBi , Aifi .— Ai -t- zBi. 

Exchanging t and s in the second integral, and replacing Aout and Aj„ by 
their expressions, an elementary computation (plus the fact that Imp = 
0(e)) leads to the required approximation of /3 i,r. 

The same procedure (but somehow simpler) shows that (^ 2,1 and (52,R 
are 0(e‘^) (note that, for and /32,r, one must use the fact that the 

integra l f^{ts)~ 2 e~^^^~^^^dtds is finite). Finally, concerning we 
see on ( |5.7| ) that it involves only functions that are real when E is real, and 
the same kind of estimates as for show that it is real up to 0(h). □ 


6. Quantization condition 


In order to simplify the writing, we will use the following notations: 
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If ui is any of the functions ^ or we set, 

ui(0) 


ui := 


cIui(O) 
0 
0 

± „,± 


and if U 2 is any of the functions or we set, 


U2 := 


0 

0 

U2{0) 

du2{0) 


With these notations, we see on Proposition 5.2 that we have, 

(6-1) _ 

Wo{E) := }V{wi^L, W2 ,l,wi^r, W2,r) = det(vi,i, vi,r, V2,l, V2,r) + 0{h), 
with. 


vi,s := + ai,5U^5 + Pi,s^t,s'^ 

V2,S := ^2,3 + “2,5U^^5 + P2,sut,s = L,R). 

Developing the determinant by multi-linearity, and observing that all the 
terms that involve at least three vectors with the same value of the index j 
vanish, we obtain, 

W2,_r) 

+ P2,LyV{Ui^R,Ui^R)y^{u2,L^'‘^2,R) 

+ 0!i^Ra2,R}V{u~^ Ul,R)'^i^2,W “ 2 , l ) 

+ ai,Ha2,L>V(u“^, uIlW{u2,w 

+ '“2,r) 

+ ai,La2,R>V(u^^, u~^) 

+ ai,La2,L>V(u+^, u+^) 

+ /5i,lVV(u^2,5 '^1,r)^(^2 ,L’ ^2 ,r) 0{h). 


In particular, we observe that each aj^s is always multiplied by another 
similar quantity, that is, by 0(s). As a consequence, an error on Oj^s of 
order will lead to a error of order h in W’o(Pl), and thus, by Proposition 
we can replace a 2 ,s by (5 = L,R). Then, computing the various 


5.3 


Wronskians that appear (see Appendix |8.3[ ), we find, 

( 6 . 2 ) 


—ITT e 


*4>Vo(ill) = - 4\/2(cos ^)(1 + C>(e^)) + 4\/2(sin 


h 


+ 8e 4 ai^Rai^L + 2\/2(sin '^)Pi,R + ^V2(sin 
-|- 2\/2(sin -j^)I3i^r + (D{h). 
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Finally, observing that we have, 
7r^ro(0)^e^ 


2 

— 


{2nAf^B + ~ ’^I^b) + C)(/i); 


ai,L«i,R = vr2ro(0)^e^\/2e*4 (sin - ifJ^B) + 0{h + e^| cos ^|); 

n n 


a{L = 47r"‘ro(0)^e^(sin + C’(/i + e"'| cos ^|); 

’ n n 

/3i,i? = Re/3i,H + i7r^ro(0)^e^(^^ + /r|) + 0{h)-, 

/3i,L = Re/?i,L + 0(/i), 

we obtain, 

(6.3) 

-z7r^e“*^>Vo(S) = - 4\/2(cos ^)(1 + 0(e^)) + 2\/2(sin ^)Re(/3i,i, + /3 i,r) 




+ 4\/2(7rro(0)e)^(sin (3/rAAiR + *(3 + sin^ w)/^A 




+ 0(/i). 


Now, by Proposition 4.2, the quantization condition reads, 

m{E) = 0. 


Hence, in view of (6.2) and of Proposition |5.2[ if we set, 

bo ■= ^Re(^gL+ +37r^ro(0)^^A(Re/9)/iR(Rep) = 0(1), 

we have proved. 

Proposition 6.1. E = pe^ G ^^{Co) is a resonance of P if and only if, 

ME) o. 


(6.4) 

where 

(6.5) 


cos 


= e^E{E,e), 


E{E, e) = 6o + ivr ro(0) (3 + sin 


2„ /n^2/'Q I „;„2 \ „^2^ ”^(-®) 


^3 -)l^A(Re/9) j sin ^3 


+ 0(e). 


Remark 6.1. The quantization condition (6.4) is of Bohr-Sommerfeld type 
associated with the single potential well ofVi{x). The imaginary part of 
E{E,e) will give an estimate on the width of resonances. 


7. Completion of the proof 


In order to solve (6.4) in Ph{Co) (where Co may actually vary a little bit in 
order to avoid ’’border” effects), we first observe that, near S = 0, the roots 
of the equation cos{A{E)/h) = 0 are given hy E = ek{h) with. 


ek{h) := A 


-1 


{k + ^)Trh 


{k G Z). 
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(Here, A: G Z must be taken such a way that A~^{{k + \)'!Th) is effectively 
close to 0.) In particular, restricting to E = 0{e‘^), and writing A{E) = 
^(0) + EA'{0) + \E‘^A!'{^) + 0{E^), we obtain the well known relation. 


ek{h) = Xk{h)e - 


2 Xk{h)‘^A^'{0) A 6 ', 


2A!{Q) 


+ 0{e^ 


where \k{h) is defined in (2.3). In particular, the distance between two 
consecutive ek{h)'s is of order h. Moreover, if i? G C is such that E = O(e^) 
and A{E) stays at a distance greater than 5h from the set {ek{h ); k G Z}, 
with d > 0 constant, then cos{A{E)/h) remains at some fix positive distance 
from 0. As a consequence, for e > 0 small enough, we can apply the Rouche 
theorem and conclude that, for each k such that Xk{h) = 0(1), there exists 
a unique solution Ei;{h) to ( |6.4[ ) such that, 

Ek(h) = ek(h) + o(h), 


and, conversely, all the roots of (6.4) in VhiCo) are of this type. 


Now, going back to equation (6.4), we immediately see that, actually, we 
have, 

Ekih) = ek{h) + Oie^), 


so that (2.4) is proved in the case ri = r 2 = 1 


In order to prove (2.5), we first observe that, by (2.4), the equation (6.4) 
implies, 

cos = e^E(Xk(h)£^,£) + O(e^). 

Then, taking the local inverse of cos near (k + \)ti, and the inverse of A 
near E = H, ( |2.5[ ) (with ti = r 2 = 1) immediately follows from (6.5) and the 
fact that sin ) = (—1)^ + 0{£^). 

When Ti and T 2 are general positive numbers, we observe that all the con¬ 
structions Sections |3.1| to and of Appendix |8.2| remain completely un¬ 
changed. Therefore, the proof proceeds exactly in the same way, and the 
only differences are in the approximate values of and ^ 2 (£~‘^y)- A 

very simple computation shows that they become, 

6(e“^y) = (^y- + O(e^); 

6(e"^y) = (y + -^^ + c>(e^). 


As a consequence, the approximations given in ( 5.10[ ) have to be changed, 
too, and indeed now the functions and fiB will depend also on the side 
where we are working (on or on J/j). On Ir (that is, for and 
they will just be as before, with Ai (y — p),Bi(y — p) substituted by 
Ai (rj“^'^^(Tiy - p)),Bi (r^'^^^iriy - p)), and Ai (y -|- p), Bi {y + p) substituted 
by Ai (t2 (T 2 y -I- p), Bi (t^ {T 2 y + p)- On II (that is, for ai^L), PA and 
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become, 

= [ Ai - p))Ai {T 2 ^^^{T 2 y + p))dy; 

J —OO 

psit) = j Bi - p))Ai (r2'^'^^(r2?/ + p))dy. 

J —OO 

But the computations proceeds in a similar way, and since PA{t) = P 2 {'t), 
the result follows in the general case, too. 


Remark 7.1. Resonant states By construction, the resonant state (fk 
associated with Ek{h) can be written both as a linear combination of wi^l 
and W 2 ,l, and a linear combination of wi^r and W 2 ,r (all computed at E = 
Ek{h)). Th e co efficients can actually be computed (up to 0{h)) by using 
Proposition 5.3 identifying each Wj^s with the vector of given by, 


^j,L 


dwJgiO) 


Then, the approximations of the various functions involved with their Airy 
representation (obtained from Propositions 8.2 and 8.5) permit us to write 
(pk as, 


(7.1) Pk = Wl,L - PW2,L = >^Wl,R + 1^W2,R, 

with, 

A = sin ^ + 0{h) = (-1)^ + O(e^); 
u = -2Aai,_R + 0{h) = 4e7rro(0)//A(AA:(/i)) + 0(h); 

/r = cos ^ - ^(/3 i,l + A/3i,ij + 1202 ,r) + 0{h) = C>(e^). 

Using the various asymptotic behaviours of the functions tCj, 5 ’s, one can 
derive the (semiclassical) asymptotic behaviour of pk in all of M. 


8. Appendix 

8.1. Inhomogeneous Airy equations. Let Ai (y) be the Airy function, 
which is characterised as the solution to the Airy equation 

(8.1) - u”(y) + yu{y) = b 

with exponential decay as y tends to + 00 , 

Ai(y) = e“i^^^''(l + 0{y~^)) [y -A- + 00 ). 

2y vr 

It is oscillating when y < 0, and its behaviour is, 

= :^(“y)~^sin^^(-y)3/2 +(1 + C>(|y|"^)) y ^ 

Moreover, the asymptotic behaviour of its derivative Ai'(y) is obtained by 
formally differentiating the previous ones, and these asymptotic behaviours 
remain valid in sufficiently small complex sectors around the real line. 
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We define another solution Bi (y) to the Airy equation by the asymptotic 
behavior as y —)■ —oo 

Bi (y) = ^ sin “ l) C>(|y|"^)) (y ^ -oo). 


Bi (y) is positive and grows exponentially for y > 0, and satisfies 


1 


-i 


Bi(y) = ^y il + 0{y )) 


From the asymptotic behaviors of Ai (y) and Bi (y) as y 
see the following properties. At first, the solutions 

(8.2) Ai (y) - iBi (y)-^ {-y)~~^ exp 


(y +oo). 

—oo, we easily 


/TT 
'KX j 4 

(8.3) Ai(y)+iBi(y) ~ (-y)“3 exp 


TT 




(y -oo); 


(y -oo). 


are outgoing and incoming respectively for negative y, and secondly, the 
wronskian of Ai (y) and Bi (y) is given by 

W[Ai, Bi] := Ai (y)Bi' (y) - Ai' (y)Bi (y) = 7r“^ 


Set 

^(y, z) := -TT {Ai (y)Bi {z) - Ai ( 2 ;)Bi (y)} , 
and define the integral operators K, and /C for / G (M) 

(8.4) K:[f]{y) -.= [ K{y,z)f{z)dz-,K[f]{y)-.= l K{-y,-z)f{z)dz. 

Jy Jo 

The function K. [/] (y) gives a particular solution to the inhomogeneous equa¬ 
tion 

-u" + yu = f, 

while K, [/] (y) gives a particular solution to 

-u -yu = f. 

Notice that there exists a symmetric property between /C and K, , namely, 

(8.5) /C[/](-y) = £[/](y), 

where /(y) = f{—y)- Moreover, if y G C, then the operators JC p and 1C p 
defined by, 

( 8 . 6 ) 

K^p[f]{y) ■= [ K{y-p,z-p)f{z)dz-, ICp[f]{y) := [ K{-y-p,-z-p)f{z)dz 
Jy Jo 

give solutions to the equations 

-u" + {y- p)u = f ; -u -{y + p)u = /, 

respectively. 

Remark 8.1. Observe that all these constructions remain valid when y 
becomes complex, as long as Im y stays bounded. In that case, the integrals 
must be taken along any complex curve joining y to 0. 
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8.2. Yafaev’s constructions. In this appendix, we recall and extend the 
constructions made in [Ya] for the scalar Schrodinger equation [Pj—E)u = 0. 
In |Yaj such constructions are made for real E only, and they just concern 
solutions decaying at infinity. Here, we need to consider complex values of 
E and exponentially large or oscillating solutions, too. 


We hx xq G (x*, 0), and we first treat the case j = 1. For E small enough, 
we denote by xi = xi{E) the only point near 0 where Vi — E vanishes (in 
particular, xi{E) = E-\-0{E'^) depends analytically on E). In the particular 
case where E is real, we can define as in |Yaj . 


^i{x;E) 

(8.7) 


^i{x;E) 



when x > xi] 


2 

3 

when xo < < ici- 


Then, it is easy to check that ^i{x] E) = x — xi{E) + 0{{x — xi)^) as 
X —)■ xi, and that depends analytically on x and Fi for x G (xo,+oo), 
E small enough. Since also Vi has a positive limit at +oo, we see that we 
can extend analytically to a complex neighbourhood of (xq, +oo) x {0}. 

Then, ^i(x) satisfies = Vi — E and Re^( > 0 everywhere on [xq, +oo). 

In particular, when E G T’/i(C'o) is fixed and x varies in (xq, +oo), then 
describes a smooth complex curve parametrised by x, with Im^i(x) = 0{h) 
uniformly. From now on, in order to simplify the notations, we drop the 
dependance of with respect to E. The result is (see also |Yaj . Theorem 
2.5), 


Proposition 8.2. Let E G T)h{Co). Then, the equation {Pi — E)u = 0 
admits two solutions ufj^ on M, such that, as x —)• + 00 , (and uniformly with 
respect to h > 0 small enough), 

u^Jx) ~ (1 + 0{h))^{Vi{x) - y/vht)-Edtlh^ 

’ vvr 

and, as /i —)■ 0+, 

^r,R(^) = 2(?i(a^))“^Ai {h-hi{x)){l + 0{h)) 

on [xo, + 00 ) n {Re|i(x) > 0}; 

u^j^ix) = 2(e)(x))-Ui {h-ki{x)) + 0{h{l + h-2/3|^i(x)|)-3)) 

on [xo, + 00 ) n {Re.^i(x) < 0}; 
= {^i{x))~^Bi{h-^i{x)){l + 0{h)) 

on [xo, + 00 ) n {Re.^i(x) > 0}; 
= (?i(^))“^Bi(/i“i^i(x)) + C>(/i(l + /i“2/^|Ci(x)|)"3)) 

on [xo, + 00 ) n {Re^i(x) < 0}. 


Proof. The proof for ^ is the same as in |Yaj . with the difference that, 
here, E may be complex. However, since we have ImFi = 0{h), all the 
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estimates in |Ya| remain valid. Observe, in particular, that when Re.^i(x) > 
0 , then Im(^i(x))^/^ = 0{h), and thus Ai{h~3^i{x)) / 0 . 

Therefore, let us focus on the construction of As in |Ya| . Section 3, 

setting t := and f{t) := |^(x) 2 u(x), the equation (Pi — E)u = 0 

becomes, 

( 8 . 8 ) - f{t) + tf{t) = Rit)f{t), 


with, 

(8.9) 


R{t) = 


pix) ■■= f(li(a:)) 


0(1 + |6(x)|)-2. 


(In the last estimate, we ha ve u sed the fact t hat \V-I{x)y + \Vy(x)\ = 0((1 + 
|x|)“^).) Defining K. as in (8.4), we reduce ( 8 . 8 ) to the Volterra equation. 


( 8 . 10 ) 


/ = Bi +/C[P/]. 


Then, a continuous solution of (8.10) will be solution of ( 8 . 8 ), too, and we 
expect it to have the right behaviour at infinity. Moreover, it is enough 
to solve (8.10) separately on Ret > 0 and Ret < 0 (where, in any case, t 
remains on the curve T := {/i“^/^^i(x); x G [xo,+oo)}). 


When Ret > 0, one has Bi (t) 7 ^ 0, and we set g := //Bi. Then g must be 
solution to. 


( 8 . 11 ) 


g = l + Lg, 


with. 


Lg{t) := TT 


Ai (t) 
Bi (t) 


Bi (s)^ — Ai (s)Bi (s) ] R{s)g{s)ds. 


Using the asymptotic behaviours of Ai and Bi at infinity, and the fact that 
f^{l + s)- ^{l + hh)-‘^ds = 0{h-i), we see that ||A||cO(rn{Ret>o}) = C’(/i). 
Therefore, (8.11) can be solved by iteration on := T n {Ret > 0}, and 
the corresponding solution to ( | 8 . 10 | ) satisfies, 

/ = Bi (t)(l + 0{h)) uniformly. 


On r_ := T n {Ret < 0}, since t = 0{h ^/^) there, and |Ai (s)| + |Bi (s)| = 
C >((|1 + |s|)“^^^) between 0 and t, we obtain 

(1 + |t|)i|/C [P/](t)| = Oih) sup (1 + \s\)-yf{s)\, 

ser_ 


and thus, ( 8 . 10 ) can be solved by iteration there, leading to a solution that 
satisfies, 


/ = Bi (t) + C>(ti(l + |t|) 4 ) uniformly. 


Moreover, the behaviour at infinity of / is obtained from that of Bi and of 
^i(x). This completes the proof of the proposition. □ 
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Remark 8.3. 

to, 


We also see on (8.11) that g' 


0{h{l + \t\Y/‘^). This leads 


=(Yi(x))'^Bi'(h ^6(a;))(l + 0(h)) 

+ 0 (/ i 2 / 3 )(/, 2/3 + (h-ki(x)) 


on [xo, +oo) n {Re^i(x) > 0}, and similar estimates are valid on [xq, +oo) H 
{Re^i(x) < 0}, and for jY)', too. For instance, on [xo,+oo) n 

{Re^i(x) < 0}, one obtains, 

(^O. L2iJ , 2 

= 2(^;(x))2Ai'(/i-3^i(x)) 

+ 0(/i5/®)(1 + (/i2/3 + |^i(x)|)-i/4). 


Remark 8.4. Similar constructions can be done on (—oo,xo], leading to 
solutions with the asymptotic behaviour. 


~ ^{Vi{x) - y/VW)-Edtlh 


where x*{E) is the only point near x* where Vi — E vanishes. 


Now, we treat the case j = 2. Here the situation is a bit different, because 
the set where V2 < 0 is unbounded, and also because there is one turning 
point only. This actually permits us to directly obtain the asymptotic of the 
solutions both at —00 and at +00. We denote by X 2 (E) the unique point 
near 0 where V 2 — E vanishes, and, when E is real, we set. 


(8.13) 


Ux;E)) 


^2(x;E) 



when X > X2; 


when X < X2. 


As before, we extend analytically this definition to complex values of E, and 
we have. 


Proposition 8.5. Let E G T>h{Co)- Then, there exist two solutions uY^ to 
equation {P 2 — E)u = 0 on M, and two constants = 1 + 0{h), such that, 
(8.14) 


u: 


± 

2,L 


(x) ~ {l+0{h))^{V2{x)-E)-^/^e'^^^FE) Edt/h 

V'TT 


(8T5) 

-02U2_^(x) ±m^M+^(x) 

~ (1 + 0{h))^{E - V2{x))-^/^e^"VETvFLdtih ^ 
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uniformly with respect to h > 0 small enough, and, 

= 2(?2(3^))"^Ai {h-h2{x)) + 0{h{l + 

OR M n {Re ^2 ( 3 ;) > 0}; 

^2 ,l(®) = 2(^2(^))~^^i {h~3^2{x)){l + 0{h)) OR M n {Re^ 2 ( 3 ;) < 0}; 

= (^2(a;))"^Bi {h-li2{x)) + 0{h{l + /i“^/^|6(a:)l)"3)) 

ORMn{Re^2(a;) >0}; 

= (^2(3;))~^Bi {h~^^2{x)){l + 0{h)) OR M n {Re,^ 2 ( 3 ;) < 0}, 

uniformly as /i —)• 0+. 


Proof. The procedure is the same as for the previous proposition, but this 
time, setting f{h~‘^^^f^ 2 {x)) ■= the equation {P 2 — E)u = 0 

becomes, 

(8.16) f{t) + tf{t) = Rit)f{t), 

where R{t) = + h‘^^^\t\)~^). In the case of ufj^, we reduce (8.16) 

to the Volterra equation. 


(8.17) 

where K. is defined by. 


/ = Ai +lC[Rf], 


(8.18) 


^[/](t):=/ K{-t,-s)f{s)ds. 


In that case, we can follow the procedure of |Yaj on {Re^2(2:) < and 

obtain a solntion with the required asymptotics (at —00 and for /i —)• 0+). 
On the other hand, on r+ := {Re.^2(3^) > 0}, we rewrite (8.17) as, 

f{t) = Ai{t) + [ K{-t,-s)R{s)f{s)ds + )C[Rf]{t), 


(where K, is as in (8.4)), and the asymptotics at infinity of Ai and Bi give, 

/C [i?/](t) = C)(/i5) sup I/I / (1 + |t|)“J(l + s)“3(l +/i5s)“^ds 

r+ Jo 

= 0{h{l + |t|)“3))sup|/|. 

r+ 

In addition, (see also |Yaj . Formula (3.18)), one has, 

/■o ht , o . 

/ K{-t,-s)R{s)f{s)ds = 0(- -^) sup I/I 

J-00 (l + |t|)5 r+ 

= 0(h(l + |t|)“a)sup|, 

r+ 


s 2 (1 + /i3s) ^ds 


Thus, (8.17) can be solved by iteration on r+, too, and there the solntion 
satisfies, 

/(t) = Ai (t) + C)(/i(l + |t|) 4)). 

The result for follows by taking := f f, 2 {x))■ 
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In the case of U 21 we reduce (8.16) to the Volterra equation, 


(8.19) 


/ = Bi +)C[Rf]. 


First working on {Ret < 0}, the same procedure as for Proposition 8.2 
leads to a solution that satisfies / = (1 + 0(/i))Bi. Next, on {Ret > 0}, the 
result follows exactly as for but this time there is no need to rewrite 
the Volterra equation. 

The asymptotic of each solution at —00 follows from that of Ai and Bi, and 
the fact that ^ 2 ix) ~ at +00 (with c > 0 constant). 

Finally, since Ai ±iBi do not vanish on M, by setting g := //(Ai ±iBi), we 
see (e. g. as in the proof of Proposition |8.2[ ) that there exist two solutions 
f± to (8.16) on M+ satisfying /± = (1 + 0{h)){Ai ± iBi). They give rise 


to two solutions v± to (P 2 — E)v = 0 (that are conjugated of each other for 
real values of E) and by computing their Wronskians with we see that 
they are of the form v± = ^ 0 / 1 x 2 ,l i with af = 1 + 0{h), so that 

(8.15) follows, too. □ 


Remark 8.6. Here again, estimates on can also be derived from 

the construction (see Remark 8.3). 


Remark 8.7. Near infinity, all these constructions extend to a complex 
sector in x, with the same asymptotic behaviours as on the real. 


Remark 8.8. These constructions can also be adapted to a problem on 
M+ (e.g., for radial solutions of a problem in MT), with potential Vj{r) 
(j = 1,2) behaving like Cjjr°^, with cj > 0 constant and 0 < a < 2, as 
r —)• 0+. In this case, —00 is replaced by r = 0, and the decaying solu¬ 
tions at —00 become solutions vanishing at 0. Thus, in the constr uction 
of uj^ (j = 1,2), it is enough to replace Ai (for instance in (8.17)) by a 
linear combination Ai + a^-Bi, where aj G M is chosen in such a way that 
Ai (—/i“^/^^j(0)) + Oj-Bi (—/i~ ^/^gj( 0)) = 0 (here ij{r) is the corresponding 
change of variable similar to (8.13)), and to use a funda mental solution of 
the Airy equation vanishing at (e.g., in (8.18), one must replace 

—00 by —h~‘^^^^j{0)). Then the proof proceeds exactly in the same way. 


For potentials behaving like Cjjr’^ at 0 (with Cj = Cj{h) > 0{, the adaptation 
is even simpler since, in that case, /j(r) —)> —00 as r —)• 0+, so that the 
construction remains the same. 
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8.3. Formulae. 





~ Bi (y - p) 

2(?i)“^Ai {h~^i) ~ 2Ai (y - p) 
(^2)“^Bi {h~^2) ~ Bi (y + p) 
2(^2) ^Ai(/i 3^2) ~ 2Ai (y + y) 




W{u2 L: U2 R 
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